Mult if r act al Properties of the Random Resistor Network 
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We study the multifractal spectrum of the current in the two-dimensional random resistor network 
at the percolation threshold. We consider two ways of applying the voltage difference: (i) two 
parallel bars, and (ii) two points. Our numerical results suggest that in the infinite system limit, the 
probability distribution behaves for small i as P(i) ~ 1/i where i is the current. As a consequence, 
the moments of i of order q < q c — do not exist and all current of value below the most probable 
one have the fractal dimension of the backbone. The backbone can thus be described in terms of 
only (i) blobs of fractal dimension dg and (ii) high current carrying bonds of fractal dimension going 
from 1/v to ds- 

PACS numbers: 64.60.Ak, 05.45.Df 



The transport properties of the percolating cluster 
have been the subject of numerous studies |Q,0j . A partic- 
ularly interesting system is the random resistor network 
(RRN), where the bonds have a random conductance. 
The random resistor network serves as a paradigm for 
many transport properties in heterogeneous systems as 
well as being a simplified model for fracture [|3|. 

The first studies of the RRN were devoted to effec- 
tive properties of the network (conductivity, permittivity, 
etc.) but for many practical applications — such as 

fracture, and dielectric breakdown [|| — the central quan- 
tity is the probability distribution P(i) of currents i. For 
instance, in the random fuse network, it is the maximum 
current corresponding the hottest or "red" bonds which 
will determine the macroscopic failure of the system || . 

The probability distribution P(i) has many interesting 
features, one of which is multifractality in order 

to describe P(i), an infinite set of exponents is needed. 
This idea of multifractality was initially proposed to treat 
turbulence fl5|| and later applied successfully in many dif- 
ferent fields, ranging from model systems such as DLA 
fL6f to physiological data such as heartbeat JTjj . 

It was first believed [|||l0| that the low current part 
of P(i) and of the multifractal spectrum follow a log- 
normal law as it is the case on hierarchical lattices. It is 
now clear fll9|] , that for small currents, the current proba- 
bility distribution follows a power law P(i) ~ where 
b > 0. For large currents, there is a weak dependence on 
the system size L. This is in contrast with small currents 
which are governed by very long paths, and therefore de- 
pend more strongly on L. It was suggested [[l8 19 that 
the exponent b of the low-current part has a 1/ log L de- 
pendence, where L is the system size. The asymptotic 
value boo of the exponent b is of crucial importance. If 
boo is finite and positive, then a low current evolves on a 
subset with a fractal dimension depending on its value. 



On the other hand, if is zero, then the low current 
part of the multifractal spectrum is flat and the entire 
backbone is contributing to low currents. It is thus im- 
portant to understand if the apparent subset structure 
with different fractal dimensions is a finite-size effect. 

This problem was adressed by Batrouni et al |T§] who 
conjectured a zero asymptotic slope and by Aharony et 
al pg j who proposed a finite asymptotic value. The max- 
imum value of L in the literature is 128 []l8|| , so numerical 
estimates could not lead to a definite conclusion. In this 
Letter, we present evidence that the asymptotic slope is 
zero. 

We first recall the basis of multifractality applied to 
the percolating two-dimensional resistor network of lin- 
ear size L. Let n(i,L) be the number of bonds carrying 
current i. By the steepest descent method, the main con- 
tribution to n(i, L) for large L is given by |]|^|l(j 



n{i,L)~LK a ^ 



(1) 



where a = — log i/ log i. The multifractal spectrum 
/(a, L) = log nj log L can thus be interpreted as the frac- 
tal dimension of the subset of bonds carrying the cur- 
rent i. The g-th moment of the current is defined as 
M q = (J2 i q ): where the sum is over all bonds carrying a 
non-zero current and (•) denotes an average over different 
disorder configurations. These moments exists for q > q Cl 
and it can be easily shown [|l9| that the "threshold" is 
q c = —b. The asymptotic slope thus give the asymptotic 
value of the threshold q c . 

For the fixed current ensemble, one observes that 



lfj M q ~ L Tq for large L and for q > q c and where T q 
is a universal exponent. In particular, To = g?_b, T2 =t/v, 
and Too = 1/f |22j where ds is the fractal dimension of 
the backbone, t the conductivity exponent, and v is the 
correlation length exponent. If the behavior is monofrac- 
tal, then r q is a linear function of q, while in the multi- 
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fractal case, the exponents are not described by a simple 
linear function of q. In the L — > oo limit, knowing /(a) 
is equivalent to knowing the infinite set of exponents r q , 
as f(a) is the Legendre transform of r q [||. 

The low current part of f(a, L) was found numerically 
to be a power law of slope b — b(L), where []18[ 

b(L) = 6oo + rr^-r + e(L) (2) 
log L 

and e(L) is a correction decreasing faster than 1/logL 
when L is increasing. This equation shows a strong finite- 
size effect since log L grows very slowly, and two possi- 
bilities for &oo were proposed, 600 = or 600 = 1/4 
& 

We consider the two-dimensional random resistor net- 
work at criticality, i.e. the fraction of conducting bonds 
p is equal to its critical value p = p c = 1/2. We first 
apply a voltage difference between two parallel bars. 
We compute f(a,L), for a fixed voltage difference, for 
L = 50, ... , 1000, and average over 10 4 configurations 
for each L. We show our results in Fig. la. The slope 
is clearly decreasing with L, confirming the strong finite 
size effects already observed ■ 

Next, we consider a second type of configuration, which 
we call the "two injection points" case, in contrast with 
the usual "parallel bars" case. We impose a voltage dif- 
ference between two points P and Q separated by a dis- 
tance r, and we look for the backbone connecting these 
two points. This situation was studied in pl| , p8[ , but 
here we keep only the backbones of size L. In this way, 
we have large backbones connecting the two points P and 
Q, and for r <C L we expect to have a large number of 
small currents on bonds belonging to long loops. The 
multifractal spectrum is then defined in the same way as 
for the parallel bars and we calculate for different values 
of L the slope of its small-current part. The multifractal 
spectrum in this case is shown in Fig. lb. We observe 
that there is a large amount of small currents, and that 
the asymptotic limit is reached faster in the two injection 
points case. We expect that the low current distribution 
will be asymptotically the same as in the parallel bar 
case, so the consistency between the two configurations 
will support our results. However, for large currents there 
are some distinct differences in the multifractal spectrum 

II- 

Fig. 2a shows the slope b versus 1/logL according to 
Eq. (||) for both multifractal spectra. The extrapolation 
to L = 00 is consistent with = in both cases. This 
result is consistent with the behavior of the successive 
intercepts (Fig. 2b). 

Another functional form of b versus L could lead to an- 
other value of boo. If we replace the abscissa of Fig. 2(a) 
by l/(logL) K , then we find that the extrapolated value 
for boo depends on k, ranging from boo ~ 0.10 for k = 2 
to boo < (which is impossible) for n = 0.5. It is nu- 
merically difficult to distinguish between a 1/ log L and 
a l/(logL) 2 behavior, but the 1/logL is the most com- 
monly used @[n§. 



In Fig. 2a, we observe higher order corrections to the 
behavior b(L) = 600 + A/logL. A better fit can be ob- 
tained by adding to the linear form a small quadratic 
term L/(logL) 2 (and eventually even cubic and quar- 
tic terms). We find that we cannot do a quadratic fit 
over the whole range of 1/logL, and indeed this leads 
to two non-physical results: (a) For both geometries, the 
fits have negative slopes at 1/logL = 0, which is not 
physical since the larger the system size, the larger the 
number of small currents, so the behavior of b should be 
monotonically decreasing with L. (b) A second defect of 
these quadratic fits is that the obtained values for the 
intercepts are different for the two geometries, which is 
impossible. 

The important assumption here is the behavior of the 
leading term. There is no proof that the leading term 
of the expansion is 1/logL rather than (l/logL) K with 
k ^ 1. However, the assumption that the leading term of 
the expansion is 1/ log L with boo = is consistent with 
our numerical data, and shows that the correction e(L) 
decays faster than an inverse power of log L (see Fig. 2c) . 

Finally, we note that the sequence of maximum values 
of f(a,L) for the two injection points case plausibly ex- 
trapolates in the variable 1 / log L as L — > 00 to a value 
of (1b close to the known value 1.64 (Fig. 3). 

Thus our results suggest the intriguing possibility that 
for L — > 00, the small current part of f(a,L) is a hori- 
zontal line at the value ds, implying that in an infinite 
system the fractal dimension of the subset contributing 
to small current is ds, independently of the value of a. 
In this sense, the small current probability distribution 
is apparently not multifractal. The "perfectly balanced" 
bonds which carry zero current have a fractal dimension 
equal to ds fl3|- Since these bonds contribute to f(a, L) 
for a — ► 00, the fact that their fractal dimension is 
supports our hypothesis that boo — 0. A related conclu- 
sion is that q c — 0, or the negative moments of the cur- 
rent do not exist in the infinite-size limit. In particular, 
it shows that the first-passage time for a tracer particle 
travelling in a flow field in a porous medium modelled by 
a percolation cluster diverges in an infinite system. 

Moreover, the result boo = is supported by the 
following argument. If boo were not zero, then the 
number of bonds carrying a small current i would be 
n{i — > 0,L = 00) ~ i b °° . This behavior would indi- 
cate that the number of bonds carrying a small current 
i approaches zero when i — > 0, which seems unlikely, 
since on an infinite backbone, the number of loops is very 
large, and n(i —> 0, L = 00) should be nonzero. Hence 
froo = 0. This argument is consistent with the fact that 
the total number of bonds carrying a nonzero current, 
L n(i, L)d{\ogi), should diverge as L — > 00. 

For large values of the current, the multifractal fea- 
tures do not change as L increase, suggesting that in the 
infinite-size limit, there are essentially two different type 
of subsets. The first comprises the blobs of fractal di- 
mension ds, and the second set comprises links carrying 
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larger values of the current (red bonds) , of fractal dimen- 
sion ranging from d rc d = 1/v to ds- 
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FIG. 1. Multifractal spectra for fixed voltage for (a) par- 
allel bars and (b) for two injection points separated by a 
distance r = 4. We show the results for 7 different values 
of L = 50, 100, 200, 400, 600, 750, 1000 and averaged over 10 4 
configurations (L increases from bottom to the top). The 
horizontal dashed line is at ds — 1.643. 
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FIG. 2. (a) Slope 6 versus 1/logL. The circles corre- 
spond to the parallel bars case and the triangles to the two 
injection points case. These values were obtained by fitting 
the small current parts of Figs. la,b, roughly over the range 
2.5 < a < 5. The extrapolation shown as a guide to the 
eye is consistent with boo = 0. The error bars where esti- 
mated by computing the local slopes and are going from 0.02 
to 0.005 as L increases, (b) Successive intercepts computed 
by using a least square fit over three successive points. The 
circles correspond to the parallel bars case, the triangles to 
the two injection points case. These plots are consistent with 
boo = 0. (c) Correction e(L) for L = 25 to L = 1000 as given 
by Eq. (2). This plot shows the fast decay of the correction. 
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FIG. 3. The maximum of f(a,L), f ma x(L), in the the 
two injection points case vs. 1/logL. The least square fit 
shown gives the extrapolated value at L — > oo of cLb close 
to the accepted value of 1.64. Also, we found that a log-log 
plot of ds — fmax(L) vs. L is remarkably straight, with slope 
~ -0.20. 
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